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Abstract 

In nuclear reactions induced by hadrons and ions of high energies, nuclei can dis- 
integrate into many fragments during a short time (~100 fm/c). This phenomenon 
known as nuclear multifragmentation was under intensive investigation last 20 years. 
It was established that multifragmentation is an universal process taking place in 
all reactions when the excitation energy transferred to nuclei is high enough, more 
than 3 MeV per nucleon, independently on the initial dynamical stage of the reac- 
tions. Very known compound nucleus decay processes (sequential evaporation and 
fission), which are usual for low energies, disappear and multifragmentation domi- 
nates at high excitation energy. For this reason, calculation of multifragmentation 
must be carried on in all cases when production of highly excited nuclei is expected, 
including spallation reactions. From the other hand, one can consider multifrag- 
mentation as manifestation of the liquid-gas phase transition in finite nuclei. This 
gives way for studying nuclear matter at subnuclear densities and for applications 
of properties of nuclear matter extracted from multifragmentation reactions in as- 
trophysics. In this contribution, the Statistical Multifragmentation Model (SMM), 
which combines the compound nucleus processes at low energies and multifragmen- 
tation at high energies, is described. The most important ingredients of the model 
are discussed. 



1 Introduction 

Statistical approaches have proved to be very successful for description of nuclear re- 
actions. According to the statistical hypothesis, initial dynamical interactions between 
nucleons lead to re-distribution of the available energy among many degrees of freedom, 
and the nuclear system evolves towards equilibrium. The most famous example of such an 
equilibrated nuclear source is the 'compound nucleus' introduced by Niels Bohr in 1936 
PP. It was clearly seen in low-energy nuclear reactions leading to excitation energies of 
a few tens of MeV. It is remarkable that the statistical concept works also for nuclear 
reactions induced by particles and ions of intermediate and high energies, when nuclei 
break-up into many fragments (multifragmentation) [2]. In the most general considera- 
tion the process may be subdivided into several stages: (1) a dynamical stage leading to 
formation of equilibrated nuclear system, (2) disassembly of the system into individual 
primary fragments, (3) de-excitation of hot primary fragments. 
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2 Formation of thermalized nuclear system 



At present, a number of dynamical models is used for description of nuclear reactions at 
intermediate energies. The Intranuclear Cascade Model was the first one used for realistic 
calculations of ensembles of highly excited residual nuclei which undergo multifragmen- 
tation, see e.g. [SI IH E]- Other more sophisticated models were also used for dynamical 
simulations of heavy- ion reactions, such as quantum molecular dynamics (QMD), Boltz- 
mann (Vlasov)-Uehling-Uhlenbeck (BUU, VUU) and other similar models (see e.g. refs. 
[6]). All dynamical models agree that the character of the dynamical evolution changes 
after a few rescatterings of incident nucleons, when high energy particles ('participants') 
leave the system. This can be seen from distributions of nucleon velocities and density 
profiles in remaining spectators [3 El E]- However, the time needed for equilibration 
and transition to the statistical description is still under debate. This time is estimated 
around or less than 100 fm/c for spectator matter, however, it slightly varies in different 
models. Parameters of the predicted equilibrated sources, i.e. their excitation energies, 
mass numbers and charges vary significantly with this time. In this case a reasonable 
strategy is to use results of the dynamical simulations as a qualitative guide line, but 
extract parameters of thermalized sources from the analysis of experimental data. In this 
case, one can avoid uncertainties of dynamical models in description of the thermalization 
process. 



3 Break-up of nuclear system into hot primary frag- 
ments 

3.1 Evolution from sequential decay to simultaneous break-up 

After dynamical formation of a thermalized source, its further evolution depends crucially 
on the excitation energy and mass number. The standard compound nucleus picture is 
valid only at low excitation energies when sequential evaporation of light particles and 
fission are the dominant decay channels [5j. Some modifications of the evaporation/fission 
approach were proposed in order to include emission of fragments heavier than a-particles, 
see e.g. [TQl EH EG]- However, the concept of the compound nucleus cannot be applied 
at high excitation energies, E* > 3 MeV/nucleon. The reason is that the time intervals 
between subsequent fragment emissions, estimated both within the evaporation models 
[T5] and from experimental data [H], become very short, of order of a few tens of fm/c. 
In this case there will be not enough time for the residual nucleus to reach equilibrium 
between subsequent emissions. Moreover, the produced fragments will be in the vicinity 
of each other and, therefore, should interact strongly. The rates of the particle emission 
calculated as for an isolated compound nucleus will not be reliable in this situation. 
There are many other theoretical arguments in favour of a simultaneous break-up at high 
excitation energy. For example, the Hartree-Fock and Thomas- Fermi calculations predict 
that the compound nucleus will be unstable at high temperatures [15]. Sophisticated 
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dynamical calculations have also shown that a nearly simultaneous break-up into many 
fragments is the only possible way for the evolution of highly-excited systems [TU] . 

On the other hand, the picture of a nearly simultaneous break-up in some freeze-out 
volume is more justified in this case. Indeed, the time scales of less than 100 fm/c are 
extracted for multifragmentation reactions from experimental data [TTl [T8]. There exist 
several analyses of experimental data, which also reject the binary decay mechanism of 
fragment production via sequential evaporation from a compound nucleus at high excita- 
tion energy. For example, this follows from the fact that a popular sequential GEMINI 
code cannot describe the multifragmentation data [T9 l l20l [2Tj . We believe that a formal 
reason of this failure is that the evaporation approaches always predict larger probabil- 
ities for emission of light particles (in particular, neutrons) than for intermediate mass 
fragments (IMFs). The GEMINI model, which assumes independent evaporation of frag- 
ments, fails also to describe angular correlations of the produced IMFs [22j. The reason 
is that at multifragmentation these correlations reflect Coulomb interaction of many frag- 
ments, but not a two-body kinematic. 

3.2 Statistical multifragmentation model 

Several versions of the statistical approach have been proposed for the description of 
multifragmentation reactions (see e.g. (2j [23J, El])- As the main de-excitation code we 
take the Statistical Multifragmentation Model (SMM), fully described in a review [2]. 
The reason is that this model was primary constructed for using after initial dynamical 
stage, and adjusted for this kind of hybrid calculations. 

The model assumes statistical equilibrium of excited nuclear system with mass number 
A , charge Z , and excitation energy (above the ground state) E at a low-density freeze- 
out volume. This volume can be parameterized as V = Vq + Vf, so the baryon density is 
p = Aq/V. Vq is the volume of the system at the normal nuclear density po ~ 0.15 fm -3 . 
Vf is the so-called free volume available for translational motion of fragments. In the 
excluded volume approximation Vf may be taken as a constant for all break-up channels, 
however, under more realistic assumption, it depends on fragment multiplicity M in the 
channels [2]. 

The model considers all break-up channels (ensemble of partitions {p}) composed of 
nucleons and excited fragments taking into account the conservation of baryon number, 
electric charge and energy. An important advantage of the SMM is that besides these 
break-up channels it includes also the compound nucleus channel, and takes into account 
competition between all channels. In this way the SMM includes the conventional evapo- 
ration and fission processes at low excitation energy, and provides natural generalization 
of the de-excitation process for high excitation energy. 

In the model light nuclei with mass number A < 4 and charge Z < 2 are treated as 
elementary stable particles with masses and spins taken from the nuclear tables ("nuclear 
gas"). Only translational degrees of freedom of these particles contribute to the entropy 
of the system. Fragments with A > 4 are treated as heated nuclear liquid drops. In this 
way one may study the nuclear liquid-gas coexistence in the freeze-out volume. Their 
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individual free energies Faz are parameterized as a sum of the bulk, surface, Coulomb 
and symmetry energy contributions 

Faz = Faz + f az + e az + e az ■ W 

The standard expressions for these terms are: F^ z = (—Wo — T 2 /e )A, where T is the 
temperature, the parameter e is related to the level density, and W = 16 MeV is the 
binding energy of infinite nuclear matter; Ff z = B$A 2 I 3 ( ^ c 2 ) 5//4 , where B = 18 MeV 
is the surface coefficient, and T c = 18 MeV is the critical temperature of infinite nuclear 
matter; E% z = cZ 2 /A 1 ^ 3 , where c = (3/5)(e 2 /r )(l — (p/po) 1//3 ) is the Coulomb parameter 
(obtained in the Wigner-Seitz approximation), with the charge unit e and r =1.17 fm; 
•^Az™ = 7^4 ~~ 2Z) 2 /y4, where 7 = 25 MeV is the symmetry energy parameter. These 
parameters are those of the Bethe-Weizsacker formula and correspond to the assumption 
of isolated fragments with normal density in the freeze-out configuration, an assumption 
found to be quite successful in many applications. It is to be expected, however, that in 
a more realistic treatment primary fragments will have to be considered not only excited 
but also expanded and still subject to a residual nuclear interaction between them. These 
effects can be accounted for in the fragment free energies by changing the corresponding 
liquid-drop parameters. The Coulomb interaction of fragments in the freeze-out volume 
is described within the Wigner-Seitz approximation (see ref. [2] for details). 

As is well known, the number of partitions of medium and heavy systems (A ~ 100) 
is enormous (see e.g. I n order to take them into account the model uses few pre- 

scriptions. At small excitation energies the standard SMM code [2] uses a microcanonical 
treatment, however, taking into account a limited number of disintegration channels: as 
a rule, only partitions with total fragment multiplicity M < 3 are considered. This is 
a very reasonable approximation at low temperature, when the compound nucleus and 
low-multiplicity channels dominate. Recently, a full microcanonical version of the SMM 
using the Markov Chain method was introduced [2H1 EE]- It can be used for exploring all 
partitions without limitation. However, it is a more time consuming approach, and it is 
used in special cases only [26J. 

Within the microcanonical ensemble the statistical weight of a partition p is calculated 

as 

W p (x exp S p , (2) 

where S p is the corresponding entropy, which depends on fragments in this partition, as 
well as on the excitation energy E , mass number A , charge Z , volume V of the system. 
In the standard treatment we follow a description which corresponds to approximate 
microcanonical ensemble. Namely, we introduce a temperature T p characterising all final 
states in each partition p. It is determined from the energy balance equation taking into 
account the total excitation energy Eq [2] . In the following we determine S p for the found 
T p by using conventional thermodynamical relations. In the standard case, it can be 
written as 

S P = ln([[ 9 a,z) + ln(H A 3 / 2 ) - ln{Al /2 ) - ln{\{ n A , z \) + 
a,z a,z a.z 
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(M - l)ln(V f /X%) + 1.5(M - 1) + - ^|^% 



where n A) z a s the number of fragments with mass A and charge Z in the partition, 
9a,z — {^sa,z + 1) is the spin degeneracy factor, X Tp = (2nh 2 / 'm N T p ) 1 ^ 2 is the nucleon 
thermal wavelength (rriN ~ 939 MeV is the average nucleon mass), and the summation is 
performed over all fragments of the partition p. We enumerate all considered partitions 
and select one of them according to its statistical weight by the Monte-Carlo method. 

At high excitation energy the standard SMM code makes a transition to the grand- 
canonical ensemble [2J, since the number of partitions with high probability becomes too 
large. In the grand canonical formulation, after integrating out translational degrees of 
freedom, one can write the mean multiplicity of nuclear fragments with A and Z as 



(n A ,z) = g A ,z^-A 3/2 exp 

Ann 



±(F AZ (T,V)-nA-vZ) 



(3) 



Here the temperature T can be found from the total energy balance of the system by 
taking into account all possible fragments with A from 1 to A and with Z from to Z 
[2J. The chemical potentials /x and v are found from the mass and charge constraints: 

J2(nA,z)A = A , ^( nAtZ )Z = Z . (4) 

A,Z A,Z 

In this case the grand canonical occupations {n A ,z) are used for Monte-Carlo sampling 
of the fragment partitions |2J. These two methods of partition generation are carefully 
adjusted to provide a smooth transition from the low energy to the high energy regimes. 



4 Propagation and de-excitation of hot fragments 

After the Monte-Carlo generation of a partition the temperature of the hot fragments, 
their excitation energy and momenta can be found from the energy balance. In this 
approach the temperature may slightly fluctuate from partition to partition, since the 
total energy energy of the system Eq is always conserved. At the next stage Coulomb 
acceleration and propagation of fragments must be taken into account. For this purpose 
the fragments are placed randomly in the freeze-out volume V (without overlapping), 
and their positions are adjusted by taking into account that their Coulomb interaction 
energy must be equal to the value calculated in the Wigner-Seitz approximation. We 
note that in the case of the Markov Chain SMM version [26J this adjustment is not 
necessary, since positions of fragments are sampled directly. In the freeze-out volume a 
possible collective flow of fragments can be also taken into account [2J . Usually it is done 
by adding additional radial velocities to the fragments (proportional to their distances 
from the centre of mass) in the beginning of Coulomb acceleration. In the following we 
resolve the Hamilton equations for motion of fragment from these initial positions in their 
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mutual Coulomb field. The energy and momentum balances are strictly respected during 
this dynamical propagation. 

The secondary de-excitation of primary hot fragments includes several mechanisms. 
For light primary fragments (with A < 16) produced in multifragmentation even a rela- 
tively small excitation energy may be comparable with their total binding energy. In this 
case we assume that the principal mechanism of de-excitation is the explosive decay of the 
excited nucleus into several smaller clusters (the Fermi break-up) [TTj 12]. In this decay 
the statistical weight of the channel p containing n particles with masses (i — 1, • • • , n) 
in volume V p can be calculated in microcanonical approximation : 
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where m = YHi=i m i * s the mass of the decaying nucleus, S = Yl^ii^Si + 1) is the 
degeneracy factor (sj is the i-th particle spin), G = Tlj=i n r * s ^ ne particle identity factor 
(rij is the number of particles of kind j). E^ in is the total kinetic energy of particles 
at infinity which can be found through the energy balance by taking into account the 
fragment excitation energy, Up is the Coulomb barrier for this decay. We have slightly 
modified this model [TT] by including fragment excited states stable with respect to the 
nucleon emission as well as some long-lived unstable nuclei. 

The successive particle emission from hot primary fragments with A > 16 is assumed 
to be their basic de-excitation mechanism, as in the case of the compound nucleus decay. 
Due to the high excitation energy of these fragments, the standard Weisskopf evaporation 
scheme was modified to take into account the heavier ejectiles up to 18 0, besides light 
particles (nucleons, d, t, a), in ground and particle-stable excited states [TT]. The width 
for the emission of a particle j from the compound nucleus (A, Z) is given by: 

Here the sum is taken over the ground and all particle-stable excited states (i = 

0, 1, • ■ ■ ,n) of the fragment j, = (2Sj + 1) is the spin degeneracy factor of the i-th 
excited state, fij and Bj are corresponding reduced mass and separation energy, E* AZ is 
the excitation energy of the initial nucleus, E is the kinetic energy of an emitted particle 
in the centre-of-mass frame. In eq. (jEJ) paz and p A > z ' are the level densities of the 
initial (A, Z) and final (A , Z ) compound nuclei. The cross section <Jj(E) of the inverse 
reaction (A , Z ) + j = (A, Z) was calculated using the optical model with nucleus-nucleus 
potential [TT]. The evaporation process was simulated by the Monte Carlo method and 
the conservation of energy and momentum was strictly controlled in each emission step. 

An important channel of de-excitation of heavy nuclei (A > 100) is fission. This 
process competes with particle emission, and it is also simulated with the Monte-Carlo 
method. Following the Bohr- Wheeler statistical approach we assume that the partial 
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width for the compound nucleus fission is proportional to the level density at the saddle 
point p S p(E) [2] : 



r '=W!-y/ Pepi^AZ — Bj — E)dE, (7) 

where Bf is the height of the fission barrier which is determined by the Myers-Swiatecki 
prescription. For approximation of p sp we used the results of the extensive analysis of 
nuclear fissility and T n /Tf branching ratios, see ref. [2] for details. 

All these models for secondary de-excitation were tested by numerical comparisons 
with experimental data on decay of compound nuclei with excitation energies less than 
2-3 MeV per nucleons. It is important that after all stages the SMM provides event by 
event simulation of the whole break-up process and allows for direct comparison with 
experimental events. 



5 Verification and applications of the SMM 

As was shown already in first publications [21 0] the SMM gives very good description 
of experimental data in the case when fragments are emitted from equilibrated sources. 
Later on, many experimental groups have successfully applied SMM for interpretation of 
their data: ALADIN (GSI, Germany) [2?], EOS (Purdue University USA) [2E], ISIS (In- 
diana University, USA) [22], Miniball-Multics (MSU, USA and INFN, Italy) [29], INDRA 
(GANIL, France) [30], FAZA (Dubna, Russia) [31], and others. In particular, the SMM 
describes charge (mass) distributions of produced fragments and their evolution with ex- 
citation energy, isotope distributions, multiplicities of produced particles and fragments, 
charge distributions of first, second, third fragments in the system, correlation functions 
(charge, angle, velocity ones) of the fragments, fragment kinetic energy distributions. Si- 
multaneously, this model reproduces global characteristics of the systems, such as caloric 
curves, critical indexes for the phase transition, different moments of the fragment charge 
distribution. Otherwords, the model can describe nearly completely experimental events, 
and, in some cases, it can be used as event generator. Importance of multifragmenta- 
tion channels for nuclear reactions is now widely recognized. The SMM is included in 
many complex codes designed to describe transport of particles and isotope production 
in matter, for example, in GEANT4 (CERN) [32]. 

From the other hand, systematic studies of multifragmentation have brought impor- 
tant information about the nuclear liquid-gas phase transition [331 [M] . Multifragmenta- 
tion reaction opens the unique possibility for investigating the phase diagram of nuclear 
matter at temperatures T « 3 — 8 MeV and densities around p ~ 0.1 — 0.3po ; which 
are expected in the freeze-out volume. Previously, only theoretical calculations without 
experimental verification were available for this phase diagram region. This information 
is crucial, for example, for construction of a reliable equation of state of stellar matter 
and modelling nuclear composition in supernovae explosions, where the same thermody- 
namical conditions of nuclear matter exist [33] . 
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One of the promising application of multifragmentation reactions is investigation of 
properties of excited nuclei embedded in surrounding of other nuclear species. They 
can be modified in comparison with the properties of isolated nuclei, since in the freeze- 
out volume fragments can interact with each other with Coulomb and residual nuclear 
forces. This study can not be performed in conventional nuclear reactions involving only 
a compound nucleus channel. However, the new properties of fragments can be extracted 
from analysis of experimental multifragmentation data. As was found there are essential 
modifications of symmetry and surface energies of hot fragments [3D1 E3 EHl EH]- These 
modified properties of fragments should be also taken into account during their secondary 
de-excitation [3D] . 

6 Conclusion 

Nearly 60 years ago nuclear physicers started to investigate the spallation reaction. Pri- 
mary it was considered as emission of few nucleons and light charged particles (and, 
eventually, fission) from a heavy nucleus. Only these processes were clearly observed that 
time, since the accelerators could provide projectile beams with relatively low energy (few 
hundred MeV). During last 15 years we have obtained solid evidences that at high projec- 
tile energies, and in heavy-ion collisions, a heavy nucleus can be completely disintegrated 
into light and intermediate mass fragments. This multifragmentation reaction is univer- 
sal, and it is a natural fast decay taking place at high energies. The multifragmentation 
channels take as much as 10-15% of the total cross section in high-energy hadron-nucleus 
reactions, and about twice more in high-energy nucleus-nucleus collisions. Moreover, 
multifragmentation reactions are responsible for production of most intermediate mass 
fragments and some specific isotopes. 

The traditional evaporation and fission models can not describe correctly this fast 
multifragmentation, since they are based on the hypothesis of a long-lived compound nu- 
cleus. There is a statistical approach, realized in the Statistical Multifragmentation Model 
(SMM), which allows for natural extension of conventional cascade-evaporation calcula- 
tions for the multifragmentation reaction. At low excitation energy it includes compound 
nucleus processes, however, at high excitations it describes the simultaneous break-up. 
Already at present the SMM demonstrates very good description of experimental data, 
especially at high excitation energy of nuclear systems (more than 3 MeV/nucleon). The 
problems, which are necessary to resolve within this approach, concern mainly a better 
description of transition from the compound nucleus to the multifragmentation decay. 
This is important for calculation of reactions initiated by low energy projectiles (up to 
1 GeV), when very few equilibrated nuclei have a high excitation energy sufficient for 
mult if r agment at ion . 
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